Final Technical Report 
NAG-1-1966 and NAG-1-1719 


Below is the final report for the aforementioned grant sponsored by your office. 

(1) Tasks Completed and Results 

a. Under the support of the NASA grant, we have developed a new geometric-optics model 
(GOM2) for the calculation of the single-scattering and polarization properties for arbitrarily 
oriented hexagonal ice crystals. From comparisons with the results computed by the finite- 
difference time domain (FDTD) method, we show that the novel geometric-optics can be applied 
to the computation of the extinction cross section and single-scattering albedo for ice crystals 
with size parameters along the minimum dimension as small as ~6. Overall agreement has also 
been obtained for the phase function when size parameters along the minimum dimension are 
larger than ~20. We demonstrate that the present model converges to the conventional ray- 
tracing method for large size parameters and produces single-scattering results close to those 
computed by the FDTD method for size parameters along the minimum dimension smaller than 
~20 (Yang and Liou 1996a). The FDTD method for the solution of light scattering by 
nonspherical particles has been developed for small ice crystals of hexagonal shapes, commonly 
occurring in cirrus clouds. In computing the scattered field in the far field and the absorption 
cross section, we have applied a new efficient algorithm involving the integration of the electric 
field over the volume inside the scatterer on the basis of electromagnetic principles. The 
applicability and the accuracy of the FDTD technique are validated by comparison with Lorenz- 
Mie scattering results for a number of size parameters and wavelengths. We demonstrate that 
neither the conventional geometric optics method nor the Lorenz-Mie theory can be used to 
approximate the scattering, absorption, and polarization features for hexagonal ice crystals with 
size parameters from approximately 5 to 20 (Yang and Liou 1996b). By unification of the 
GOM2 and FDTD methods for size parameters larger and smaller than about 20, we now have 
the theoretical and numerical tool to compute the single-scattering and polarization properties of 
ice crystals for all sizes and shapes that can be numerically defined, such as solid and hollow 
columns, double and single plates, bullet rosettes, dendrites, aggregates, and ice particles with 
rough surfaces. This is referred to as the unified theory for light scattering by ice crystals. 

b. On the satellite remote sensing algorithm development and validation, we have developed 
a numerical scheme to identify multimultilayer cirrus cloud systems using AVHRR data. It is 
based on the physical properties of the AVHRR channels 1-2 reflectance ratios, the brightness 
temperature differences between channels 4 and 5, and the channel 4 brightness temperatures. In 
this scheme, clear pixels are first separated from cloudy pixels, which are then classified into 
three types: cirrus, cirrus/low cloud, and low clouds. We have applied this scheme to the 
satellite data collected over the FIRE II IFO area during nine overpasses within seven 
observation dates. Determination of the threshold values used in the detection scheme are based 
on statistical analyses of these satellite data. Validation of the detection results has been 
performed against the cloudy conditions inferred from the collocated and coincident ground- 
based lidar and radar images, balloonbome replicator data, and NCAR CLASS humidity 
soundings on a case-by-case basis. In every case, the satellite detection results are consistent 
with the cloudy conditions inferred from these independent and complementary measurements. 
The present scheme is well suited for the detection of midlatitude multilayer cirrus cloud systems 
and tropical anvils. 
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3. Graduate Students who Participated: P. Yang: P. Rolland 

I hope that you find the aforementioned report appropriate and acceptable. Thank you for your 
support of my research program at UCLA. 


